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A new LDA+U band calculation is performed on the Kondo insulator material YbBi2 and 
an energy gap of about O.OOlRyd is obtained. Based on this, a simple tight-binding model 
with 5d£ and 4f Fs orbitals on Yb atoms and the nearest neighbor cr-bonds between them is 
constructed with a good agreement to the above the LDA+U calculation near the gap. The 
density of states is also calculated and the shape is found to be very asymmetric with respect to 
the gap. A formation mechanism of the gap is clarified for the first time in a realistic situation 
with the orbital degeneracies in both conduction bands and the f states. This model can be a 
useful starting point for incorporating the strong correlation effect, and for understanding all 
the thermal, thermoelectric, transport and magnetic properties of YbBi2. 
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§1. Introduction 

Among the rare-earth compounds which include the strongly correlated f electrons, SmBg,^) 
YbBi2,2) Ce3Bi4Pt3,^) CeRhSb,^) CeFe4Pi2,^^ CeNiSn,'^) TmSe^) etc. exhibit insulating behavior 
at low temperatures whereas the Kondo-like behaviors (the enhanced electric specific heat and 
the enhanced paramagnetic susceptibility, etc.) are often observed at higher temperatures. These 
materials are called Kondo insulators, Kondo semiconductors or heavy fermion semiconductors.^^ 
Recent studies on CeNiSn uncovered that this material is a semimetal with a pseudogap.^^ TmSe 
orders antiferromagnetically below Tn=5K, so that the insulating behavior seems to be due to 
the gap caused by this order. Sm has a complicated f shell with five to six f electrons (although 
the ground state of has vanishing total angular momentum, so that the treatment of Sm ions 
might not be so difficult than expected^'^^ ) . We exclude these three materials from the following 
discussions. 
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The Kondo effect usually occurs in metals, so that the existence of an energy gap at the Fermi 
energy in insulators may seem contradictory to the occurence of the Kondo effect in these materials. 
However, it is already clarified that the Kondo effect can occur if the Kondo temperature Tk exceeds 
the gap size Eg}^^ Even in the compounds with the periodic array of rare-earth ions, the effect of 
the strong correlation is mainly to renormalize the gap size if the band calculation does yield an 
energy gap.^^^ Thus, we can understand that the Kondo insulator is a band insulator with a strong 
correlation.^^' In fact, all the known Kondo insulators have even number of electrons which fill 
the bands below the gap. (TmSe has odd number of electrons, so that it may not be classified into 
the Kondo insulators also in the present sense.) In this context, the concept of the Kondo insulator 
may be extended beyond the rare-earth compounds. FeSi^^^ is considered to be such an example 
among the transition metal compounds with 3d electrons^^^ although the correlation may not be 
so strong as in rare-earth compounds. 

Recently, some of the Kondo insulator compounds are attracting renewed interests because of a 
possibility to be a potential candidate for an efficient thermoelectric device.^^^ However, quantitative 
analysis on these materials have been hindered because of the lack of a simple description of the 
basic electronic structures. LDA band calculations on the Kondo insulators including the f-electrons 
as itinerant ones are carried out, giving rising to an energy gap around the Fermi energy, or at 
least a tendency towards the opening of the gap.^^^ Nevertheless, the obtained band structures look 
rather complicated, so that a simpler tight-binding model description is necessary to explore the 
effect of the strong correlation starting from the band structure calculations. 

The simplest theoretical model to describe the Kondo insulators is the periodic Anderson model 
(PAM). Up to now, mostly the case with only the spin degeneracy has been investigated.^^' 
In this case, the conduction band dispersion ejj. and the f level depicted in Fig. 1(a) mix up to yield 
the bands shown in Fig. 1(b) after the hybridization, and the energy gap opens. However, it has 
already been criticized that this too simple scenario does not work in realistic systems with orbital 
degeneracy. It was pointed out by Anderson^^) that only one orbital among f states can mix with 
the plane wave state, and other f states remain unmixed. This situation was treated later in more 
detail.^^'^^) Next, consider the case shown in Fig.2(a) where the conduction band has the two-fold 
orbital degeneracy. (The spin degeneracy is not explicitly shown here and henceforth.) It is clear 
from Fig.2(b) that the bands after the hybridization can not have a gap.^^^ 

Sometimes, a PAM with orbital degeneracy in both conduction bands and f states is consid- 
ered.^^' In this case, however, an artificial assumption that the local symmetry of the f state 
is conserved even when an electron propagates from site to site, is made in order to facilitate the 
calculation. 

Theoretical studies have been performed also on the model in which the conduction band is 
treated by the free-electron model but the anisotropic hybridization matrix elements with the f- 
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electron states under the crystalline electric field (CEF) are taken into account. ^^"^'^^ Applicability 
of this type of model may, however, be limited to the materials which have such a free-electron-like 
conduction band. Thus, we need a more sophisticated model and an explanation of the formation 
of a gap in real materials based on reliable band calculations. 

(a) (b) 




Fig. 1. The simplest picture for the formation of the hybridization gap. The dispersion curves (a) without mixing 
and (b) after the mixing is introduced. 



(a) Et (b) Et 




Fig. 2. The case that the hybridization gap does not open because of the degeneracy in the conduction band. The 
dispersion curves (a) without mixing and (b) after the mixing is introduced. 

YbBi2 is the best material to study the formation mechanism of the gap because it is one of 
the most intensively investigated materials. Thermal, '^^^ transport, '^^^ optical, PES,^"^^ neutron 
scattering'^^-* experiments have been performed. In addition, the crystal structure is the simple 
NaCl type, where Yb ions and B12 clusters are located at the interpenetrating fee sites. Previous 
LDA band calculation, ""^^^ however, resulted in a semimetal with small overlap of the conduction 
and valence bands. Role of the spin-orbit interaction was found to be important by the comparison 
with the calculation for LuBi2 without it."'^) We will report a new calculation using the LDA+U 
method^'''^ in the next section and find that a gap does open. To clarify the mechanism of the 
opening of the gap when the orbital degeneracy exists, we will construct a simple tight-binding 
model for the conduction band and introduce a hybridization with the 4f Fg state in §3 and the 
formation of the hybridization gap is discussed. §4 is devoted to classifying and clarifying the 
formation mechanism of the gap. Conclusions will be given in §5. 
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§2. LDA+U Calculation for YbBi2 

We have performed a new FLAPW band calculation for YbBi2 using the LDA+U method^^^ 
by adding the term H' = Tjimm'aa' K"^'^)^m' (^"^'"^'1 *° LDA Hamiltonian, where = 
— nj^,), i7 is a parameter and is chosen to be 0.5 Ryd, and = X]nk(^"^<^l^^)(^^N'^''^') 

is the density matrix between the orbitals m, m! and the spins a, a' . Here, \nk) is the Bloch function, 
and we consider only ^ = 3 for the localized orbital \ima) at site i. Note that the original LDA+U 
method is extended to the case with the spin-orbit interaction by introducing the off-diagonal 
matrix elements with respect to the spin indices in nj^, which is calculated self-consistently. 

Introduction of the potential V^^, pulls down the occupied 4f levels compared to the LDA 
positions. In the case of YbBi2, the resultant 4f level lies below the 5d conduction band. We 
think that it is an artifact of the LDA scheme, in which the exchange energy may be overcounted 
for the 4f electons of heavy Yb atom compared to the 2p electrons of the light B atom, because 
the 4f electrons are located around the core-electrons which contribute the LDA potential. This 
phenomenon can not be corrected even by introducting the LDA+U method. Therefore, we have 
pulled up the 4f levels additionally by 0.3 Ryd. Whole the calculation was made self-consistent 
including this additional shift, and then we obtained the 4f levels located in the conduction band, 
which has t2g character and consists mainly of 5de states on Yb. We consider that these positions 
of the 4f levels are reasonable since otherwise we can not obtain the gap. The 4f levels are split 
into Fe, Fy and Fg under the cubic symmetry with the energies Ey^ < Ey-j < Ey^. 

We found that a small gap of about 0.0013 Ryd opens due to the mixing of the 4f Fg state 
and the conduction band of the t2g character as shown in Figs. 3 and 4. Number of f electrons is 
about 13.3. Note that F7 and Fg states are lower in energy than Fg, whereas the recent neutron 
scattering experiment^^^ found the opposite crystalline field level scheme Ey^ > Eyj > Ey^. This 
contradiction can be simply understood by taking the hole picture. In this view point, the ground 
state of 4f^^ corresponds to the one hole state in Fg, and the excitation to Fy or Fq states needs 
to annihilate this hole and create a hole in F7 or Fg states, which needs positive energy in total. 
The density of states is displayed in Fig.5. Note that the shape is very asymmetric with respect 
to the gap. Whole of the electronic band structure seems very complicated, but if one looks at the 
dispersion curves near the gap, one finds that it is not so much complicated. We will express them 
by a simple tight-binding model in the next section. 



§3. Tight-binding Model for YbBi2 

We express the LDA+U bands near the gap by the simple tight-binding model'^'^) with only the 
(ddcr) overlapping integral between 5ds {xy, yz and zx) and 5d7 (x^ — and 3z^ — r^) orbitals on 
Yb ions. This (ddcr) should be regarded as being produced by the effective hopping through the 
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Fig. 3. LDA+U band calculation for YbBi2. The three flat bands correspond to the 4f Fg, F7 and Fg states, 
respectively. 

B12 clusters. Configuration of the orbitals in (001) plane of the crystal is displayed in Fig. 6. We 
locate the energy levels of 5de and 5d7 orbitals at Ef^^=1.0 Ryd and E^^ =1.4 Ryd, respectively, 
and set {dda)= 0.06 Ryd. Usually, (ddfi) is negative, but it is set positive here since the hopping 
through B12 clusters may change the sign, and only this choice can reproduce the LDA+U band 
calculation. The resulting bands, shown in Fig. 7 along the symmetry axes r(000)-K(|||)-X(110)- 
W(1^0)-L(iii)-r(000)-X(100), consist of the lower 5de and the higher 5d7 bands, and look very 
similar to the dispersions of the conduction bands of i2g and Eg characters near Ep in the LDA and 
LDA+U calculations, if the f states are removed. 

Furthermore, if we take only the de orbitals, the Hamiltonian matrix becomes diagonal and the 
energy dispersions are given by the following simple expressions: 

Ef = E,, + 3(dda) cos(^^) cos(^), (1) 

where {a,/3) = (x,y), {y,z) and {z,x). The splitting at L point in Fig. 7 disappeared in this 
simplified model. It is easy to see that the lowest band along r-K-X(llO) in Fig. 7 is doubly 
degenerate and given by E^ and E^ with kz = 0, whereas the next band has no degeneracy 
and is given by E!^ with kx = ky and kz = 0. The lowest bands along r-X(lOO) are also doubly 
degenerated. 

It is expected that the inclusion of the f states at 0.88 Ryd and the mixing with the conduction 
bands in Fig. 7 will yield the bands with the energy gap as obtained by the LDA+U calculation. 
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Fig. 4. Expanded figure of the dispersion curves near the energy gap. 




Fig. 5. The density of states of YbBi2 calculated by the LDA+U method. The insert shows the expanded view of 
the bands near the gap. 



However, the lowest branch in Fig. 7 along F-X line is doubly degenerate, so that a simple mixing 
band picture cannot yield an energy gap, as was mentioned in §1. One has to take account of 
the symmetry and degeneracy of the 4f states with spin-orbit interaction (Fg) properly as well as 
those of the conduction bands (t2g)- The band calculation does include these features and yield 
the energy gap. Therefore, our tight-binding band will also yield an energy gap if these features 
are taken into account. 

The mixing matrix elements (Slater-Koster integrals) between d and f states have been given 



6 



Fig. 6. Configuration of the de orbitals (large) on Yb sites and p orbitals (small) on B sites in (001) plane. The 
shadows indicate the negative-value parts of the wavefunctions. 




Fig. 7. Energy dispersion of the tight-binding model for the YbBi2 conduction band. The primary orbital xy, yz 
or zx in each band is designated for the lower de bands. 



only for the d and f states under cubic CEF and without spin-orbit interaction.^^) We need an 
expression including spin-orbit interaction. The Fg states under cubic CEF in the subspace of the 
total angular momentum J = 7/2 (for Yb) are expressed in terms of the spherical harmonics Yp^s 



with £ = 3 and the spinors x±'s as. 
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Yg^'s can be related to the normalized cubic harmonics by 
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Using these relations, one obtains 



l^[±\x{5x' - 3r')) - i\y{5y' - 3r'))} 
-^y^{±|x(5x2 - 3r2)) - i|y(5y2 - Sr^))} 



(4b) 



The mixing matrix elements between the above cubic harmonics and d orbitals under cubic CEF 
are given in ref.41. After lengthy calculations, the matrix elements between d and f states at k are 
finally given by 
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where Cq = cos(/ca/2), = sin(/ca/2) (a = x,y,z) and t = (dfcr). Note that we have retained only 
the nearest neighbor (dfu) bonds as the simplest model. Diagonalizing the Hamiltonian matrix 
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for E^g = 0.88 Ryd, (dd(T)=0.06 Ryd and (df(T)=0.01 Ryd, we found that a gap almost opens 
but a small overlap of the bands remains between W and L points. It was not improved by the 
introduction of (dfTr), (ffcr) nor the second nearest bonds. Therefore, we have shifted down the 
bands below the gap by AE = —0.005 Ryd relative to the bands above the gap, similar to the 
LDA+U treatment, and obtain the dispersion curves shown in Fig.8 which have an indirect gap of 
about 0.003 Ryd between X and L points. The integrals (dfTr) = -0.005 Ryd and (ffa) = -0.002 
Ryd are also included here to improve the agreement with the band calculation. Note that (ffcr) 
lifts the orbital degeneracy of Fg states except the symmetry points. The value of the gap is larger 
than that in experiment, but we used it for demonstrating the gap clearly. We also expect that 
it will be renormalized down to a smaller value by the correlation effect, although the gap size is 
of course adjustable here. The density of states is shown in Fig.9, which is very asymmetric with 
respect to the gap, similar to that obtained by the LDA+U method. 



1.2 




Fig. 8. Energy dispersion of the tight-binding model for the YbBi2 including the hybridization with 4f Fs states. 



§4. Formation Mechanism of Hybridization Gap 

In order to understand the mechanism of the opening of the gap, we first consider a simple example 
with degeneracy in both extended and localized orbitals. Fg states have four-fold degeneracy, which 
consists of two-fold orbital degeneracy times two-fold Kramers one. The Kramers degeneracy always 
remains even after the hybridization or the inclusion of the spin-orbit interaction. Therefore, we 
take xy, yz and zx orbitals as the extended states and regard — y^ and Sz^ — as the localized 
ones to imitate the two-fold orbital degeneracy of Fg states, and construct the tight-binding model 
with nearest-neighbor (ddcr) bonds only between de-de and d£-d7. The Hamiltonian matrix reads 
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Fig. 9. Density of states of the tight-binding model for the YbBi2 band. 
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(6) 



are given by eq.(l). Usually, t^j must be equal to (ddcr) but here we assume t^^ ^ {dda). In this 
model, the mixing matrix elements are proportional to sin(/cj^/2) sin(/s2/2) and sin(A;2/2) sin(A;a;/2) 
for the lowest yz and zx bands, respectively, so that the gap vanishes along the symmetry axes 
r-K-X(110)-W(liO) and r-X(lOO) where = 0, as shown in Fig.ll. The density of states is shown 
in Fig. 12 and becomes singular at the flat band. Note that a gap does open, of course, if we replace 
sin(fca/2) by a constant value in eq.(6), but it does not open if we replace all the off-diagonal matrix 
elements between d7's and de's by the same value of a constant. This is because of an accidental 
degeneracy, and can be more easily seen by considering the simpler Hamiltonian matrix 
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(7) 



which has the eigenvalues e^, Ef and \£\^ + Ef ± ^(£k - Eff + 16F2]/2. The dispersion curves 
look like Fig. 10. 
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Fig. 10. The dispersion curves corresponding to the Hamiltonian matrix eq.(7) is shown. 



Note also that if the Hamiltonian matrix is given by 
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(8) 



the eigenvalues are the hybridized bands [£k + -E'/=t y (^k ~ Ef Y + 161^^] /2, each of which is doubly 
degenerated and has a gap. This or similar models have been often used in the literature,^^'^^^ but 
are not appropriate since the local symmetry of the f state is conserved even through the hopping 
from one site to the other. 

Based on these analyses, we consider eq.(5) again. The structure of this matrix looks rather 
similar to that of eq.(6), although the inclusion of the spin-orbit interaction makes it complex. If 
we set all Sq = Cq, = 1 in eq.(5), then we obtain a gap of about 0.003 Ryd without an LDA+U shift 
/S.E (Fig. 13). On the other hand, gap does not open because of accidental degeneracy if we give 
the same value of a constant to all the off-diagonal matrix elements. 

We can summarize the formation mechanism of an energy gap in the Kondo insulators as follows. 
As shown in Fig. 14, if there is a two- fold orbital degeneracy (besides the Kramers one) in the f states 
and no degeneracy in the crossing conduction band, a gap can open, but there is a possibility that 
the bands might overlap due to some dispersion of the localized states. When both the conduction 
band and the f state have two-fold degeneracy, a gap can also open as shown in Fig. 15 except the 
case of accidental degeneracy or the strong dispersion in f state. Furtheremore, complex off-diagonal 
matrix elements due to the spin-orbit interaction makes it difficult for the gap to close. 

On the other hand, gap does not open if the f state has no orbital degeneracy as in Fg or Fy states 
and the conduction band has degeneracy (Fig.2). Therefore, the experimental fact that there is a 
gap in YbBi2 indicates that the ground state is Fg. (Note if F7, say, is above the gap, and Fg is at 
the gap, the gap opens, but the number of f electrons in the ground state becomes twelve, which 
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Fig. 11. The tight-binding band with de-de and de-d7 integrals only. The parameters are chosen as tc-y = 0.5(ddcr). 



200 




Fig. 12. The density of states of the tight-binding band shown in Fig.ll. 



can not be accepted based on the experimental facts.) 



§5. Conclusions 

We have performed a new LDA+U band calculation for the most typical Kondo insulator YbBi2 
and obtained a gap of about 0.0013 Ryd. Based on this calculation, we constructed a simple tight- 
binding band model to express the bands near the energy gap. The conduction band consists of 
the 5d£ orbitals on Yb and the effective overlap integral (ddcr) is regarded as being produced by 
the hopping through the B12 clusters. This model can describe the t2g conduction band very well. 
Inclusion of the mixing with the 4f Fg states resulted in the energy gap after additional shift of the 



12 




Fig. 13. Energy dispersion of the tight-binding model for the YbBi2 in which all Sa's and Ca's are replaced with 
unity. This calculation does not include (dfyr) nor (ffcr). 




Fig. 14. The case that the hybridization gap does open. The dispersion curves (a) without mixing and (b) after 
the mixing is introduced. 




Fig. 15. The case that the hybridization gap does open. The dispersion curves (a) without mixing and (b) after 
the mixing is introduced. 



filled bands similar to the LDA+U treatment. The formation mechanism of the gap in the realistic 
situation with the degeneracy of the conduction bands and the f states is classified and clarified for 
the first time. 

The present model is very useful in constructing a theory with correlation effect, and may consis- 
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tently explain all the thermal, thermoelectric, transport and magnetic properties of YbBi2. Such 
a calculation is now in progress. Preliminary calculation for the thermoelectric power is already 
reported,^^) in which the asymmetry of the density of states with respect to the gap plays an im- 
portant role in the temperature dependence of the thermopower. Finally, it should be emphasized 
that the f-f hopping is also important, although small, to reproduce the band calculation, so that a 
use of a too simplified model with the plane wave conduction electrons and the completely localized 
f-electrons may need take care. Especially, anomalies in the density of states^^"^°) obtained for such 
a simplified model may be sensitive to the inclusion of the f-f hopping. 
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